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Abstract. We describe the weak localization correction to conductivity in ultra- 
thin graphene films, taking into account disorder scattering and the influence 
of trigonal warping of the Fermi surface. A possible manifestation of the chiral 
nature of electrons in the localization properties is hampered by trigonal warp- 
ing, resulting in a suppression of the weak anti-localization effect in monolayer 
graphene and of weak localization in bilayer graphene. Intervalley scattering due 
to atomically sharp scatterers in a realistic graphene sheet or by edges in a narrow 
wire tends to restore weak localization resulting in negative magnetoresistance in 
both materials. 

1 Introduction 

The chiral nature of quasiparticles in ultra-thin graphitic films [1I2I3I4I5] recently revealed 
in Shubnikov de Haas and quantum Hall effect measurements 6 7 8 9] originates from the 
hexagonal lattice structure of a monolayer of graphite (graphene) . The low energy behavior of 
monolayer graphene is explained in terms of two valleys of Dirac-like chiral quasiparticles with 
'isospin' linked to the momentum direction, exhibiting Berry phase ir |l|2|3|4j . Remarkably, the 
dominant low energy quasiparticles in a bilayer are different: massive chiral quasiparticles with 
a parabolic dispersion and Berry phase 2tt [5]. 

In existing graphene structures, scattering occurs predominantly from potential perturba- 
tions which are smooth on the scale of the lattice constant a. This smooth potential arises 
from charges located in the substrate at a distance d from the 2D sheet, a -C d < h/pp (h/pp 
being the Fermi wavelength) . Such a smooth potential is unable to change the isospin of chiral 
electrons so that, in a monolayer, there is a complete suppression of electron backscattering 
from potential disorder [4110] . In the theory of quantum transport in disordered systems [11] 
the suppression of backscattering is known as the anti-localization (WAL) effect [12] and, in 
monolayer graphene with purely potential scattering, a possible WAL behavior of conductivity 
[10113114115) has been related to the Berry phase ir specific to the Dirac-like Hamiltonian. Ow- 
ing to the different degree of chirality in bilayer graphene, related to Berry phase 2ir [5], purely 
potential scattering would have a different effect: no suppression of backscattering leading to 
conventional weak localization (WL) [14116] . 

In realistic graphene, there are other considerations that appear, at first glance, to be merely 
small perturbations to this picture, but they have a profound impact on the localization prop- 
erties if their effect is perceptible on length scales less than the phase coherence length. This in- 
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Fig. 1. (a) Typical magnetoresistance behavior expected in a phase-coherent (r<p ^> n) monolayer of 
graphene for a weak intervalley scattering, r* <C Ti (solid line) and for the case when the symmetry- 
breaking intravalley scattering is slower than the intervalley one r* 3> Ti (dashed). In both cases, we 
assume that the phase coherence time determines the longest relaxation time scale in the system, (b) 
Magnetoresistance of bilayer graphene, t* <C Tj (solid line) and r* 3> tj (dashed). Note that, for the 
case when r v <C Ti and r, -C Ti, Ap(B) = in both monolayer and bilayer graphene. 



eludes influence of ripples on the graphene sheet [17j leading to a weak randomization of carbon 
7r-bands, scattering off short-range defects that do not conserve isospin and valley |10ll3ll4ll5j . 
and trigonal warping of the electronic band structure which introduces asymmetry in the shape 
of the Fermi surface about each valley |4ll5ll6j . Both of them tend to destroy the manifesta- 
tion of chirality in the localization properties, resulting in a suppression of the WAL effect in 
monolayer graphene [T7] and of WL in bilayers [TS]. Moreover, owing to the inverted chirality 
of quasiparticles in the two valleys, intervalley scattering will wash out any Berry phase effect 
and restore conventional weak localization (WL) behavior of electrons in both monolayers and 
bilayers in the regime of long- lasting phase coherence [10113114115116117118119] . 

Two typical magnetoresistance curves for monolayer graphene are sketched in figure [lja). 
They illustrate two extremes: r* 3> Ti and r* <C T; where r; is the intervalley scattering time 
and t* is the combined scattering time of intravalley and intervalley scattering and of trigonal 
warping (see equation ([15|16| and ([27) below). When r* ^> r;, the magnetoresistance p(B)—p(0) 
changes sign at the field Bi such that tb ~ iy. from negative at B < Bi to positive at higher 
fields. This behavior resembles the low-to-high field crossover in the quantum correction to the 
conductivity of metals with strong spin-orbit coupling, though with an inverted sign of the 
effect. In the case of r* <C Tj, the magnetoresistance is typically of a WL type, with almost no 
sign of anti-localization up to the highest fields, which shows that, unlike in a ballistic regime or 
a quantizing magnetic field |3I5I20| . the chiral nature of quasiparticles does not manifest itself 
in the weak field magnetoresistance of realistic graphene structures. In bilayer, however, slight 
enhancement of WL behavior is expected in the case of weak intravalley symmetry breaking 
scattering, r* -C t- 1i due to different Berry phase 2tt. In the case of very strong intravalley 
symmetry breaking scattering, t* 3> Ti conventional WL magnetoresistance is expected. 

The WL behavior in graphene is novel because, with the exception of spin-orbit coupling 
[12|21j , qualitative features of WL do not usually depend on the detail of the electronic band 
structure and crystalline symmetry. In gapful multi- valley semiconductors only the size of WL 
effect may depend on the number of valleys and the strength of intervalley scattering [22123124] . 
The low-field MR, Ap{B) = p{B) — p(0), in a two dimensional electron gas or a thin metallic 
film [Illl2l22l25j in the absence of spin-orbit coupling is characterized by 



Here F(z) = lnz + ip(h + i), t v is the coherence time, D is the diffusion coefficient, and the 
integer factor se depends on whether or not states in n v valleys are mixed by disorder. This 
factor is controlled by the ratio 6 — Ti/r v between the intervalley scattering time Tj and the 
coherence time t v . In materials such as Mg, ZnO, Si, Ge, listed in Table I, where each of the 
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Table 1. Weak localization factor sg in conductors with a multi- valley conduction band and negligible 
spin-orbit coupling. The factor sg is specified for two limiting cases, no inter-valley scattering 6 = 
Ti/r^ — * oo, and for fast inter-valley scattering 9 — » 0. 
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Fermi surface pockets is p — > — p symmetric, intervalley scattering reduces the size of the WL 
MR from that described by Soo = 2n v when = t-Jt v 3> 1 to sq = 2 for S < 1. 

A more interesting scenario develops in a multi-valley semimetal, where the localization 
properties can be influenced by the absence of p — > — p symmetry of the electronic dispersion 
within a single valley, and graphene is an example of such a system. Here, we demonstrate 
how the asymmetry in the shape of the Fermi surface in each of its two valleys determines the 
observable WL behavior sketched in figure [Ha,b) ■ It has a tendency opposite to that known 
in usual semiconductors and metals: a complete absence of WL MR for infinite t; (soq = 0) 
and the standard WL effect in the limit of r; <C r v (sq — 2). In Section [5] we describe the WL 
effect in monolayer graphene with a description of the low energy Hamiltonian in Section 12.11 
a qualitative account of interference effects in Section l2~2l the model of disorder in Section [2~3l 
an account of our diagrammatic calculation of the weak localization correction in Section 12.41 
and the resulting magnetoresistance in Section 12.51 Section [3] describes the weak localization 
correction and magnetoresistance in bilayer graphene. 



2 Weak localization magnetoresistance in disordered monolayer graphene 
2.1 Low energy Hamiltonian of clean monolayer graphene 

The hexagonal lattice of monolayer graphene contains two non-equivalent sites A and B in the 
unit cell, as shown in figure [DJa). The Fermi level in a neutral graphene sheet is pinned near 
the corners of the hexagonal Brillouin zone with wave vectors K± = ±(|/ia _1 ,0) where a is 
the lattice constant. The Brillouin zone corners K± determine two non-equivalent valleys in 
the quasiparticle spectrum described by the Hamiltonian [1I4I15I28] . 

Hi = vH z (a x p x + a-ypy) + h lw , (2) 
h lw = ^IIq [a y (p x p y + Pyp x ) - a x (p 2 x - p 2 y )] . 

This Hamiltonian operates in the space of four-component wave functions, <P — [0k + (A), 4>k. + (-B), 
0k_ (£>), 4>K- (A)) describing electronic amplitudes on A and B sites and in the valleys K-|-. Here, 
we use a direct product of 'isospin' (AB lattice space) matrices <r = 1, & x , y ,z and 'pseudospin' 
inter/intra- valley matrices IIq = l,n XtVtZ to highlight the difference between the form of H\ 
in the non-equivalent valleys. The Hamiltonian H\ takes into account nearest neighbor A/B 
hopping in the plane with the first (second) term representing the first (second) order term in 
an expansion with respect to momentum p measured from the center of the valley K.±. 

Near the center of the valley K + , the Dirac-type part, v a ■ p, of Hi determines the linear 
dispersion e = vp for the electron in the conduction band and e = — vp for the valence band. 
Electrons in the conduction and valence band also differ by the isospin projection onto the 
direction of their momentum (chirality): a-p/p = 1 in the conduction band, a -p/p = — lin the 
valence band. In the valley K_, the electron chirality is mirror- reflected: it fixes a ■ p/p = — 1 
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Fig. 2. (a) Schematic plan view of the monolayer lattice containing two sites in the unit cell, A (white 
circles) and B (grey), arranged on an hexagonal lattice (solid lines). (b)Fermi lines (solid lines) in the 
vicinity of two inequivalent valleys K+ and K_ of the hexagonal Brillouin zone (dashed line). Trigonal 
warping produces asymmetry of the dispersion at each valley e(K±,p) 7^ e(K±, — p), where momentum 
p is determined with respect to the center of the valley, but the effects of warping in the valleys have 
opposite signs, e(K±,p) = e(K=p, — p). 

for the conduction band and a ■ p/p = 1 for the valence band. For an electron in the conduction 
band, the plane wave state is 

<Z>K ± , P = (±e-^ 2 | T)k ± , p + e^ 2 \ |) K± ,p) , (3) 

*K ±) _ P = (Te-^ 2 ] T> K± ,-p + e*l*\ |) K± ,-p) • (4) 

Here | T)k + , p = [1,0,0,0], | |) K+ , P - [0,1,0,0] and | T)k_, p = [0,0,1,0], | T)k_, p = [0,0,0,1], 
and the factors e ±lip ^ 2 take into account the chirality, with angle ip defining the direction of 
momentum in the plane p = (p cos ip , p sin tp) . The angular dependence w(p) ~ cos 2 (<^>/2) of the 
scattering probability off a short range potential which conserves isospin is shown in figure [3Ja). 
It demonstrates the fact that the chiral states Eqs. (|3I4) with isospin fixed to the direction of 
momentum display an absence of back scattering [4 29; 10], leading to a transport time longer 
than the scattering time r^ r = 2tq. 

The term hi-m in equation can be treated as a perturbation leading to a trigonal deforma- 
tion of a single-connected Fermi line and p — > — p asymmetry of the electron dispersion inside 
each valley illustrated in figure[2^b): e(K±, p) ^ e(K±, — p). However, due to time-reversal sym- 
metry [30] trigonal warping has opposite signs in the two valleys and e(K±,p) = e(K T , — p). 
The interplay between the two terms in H± resulting in the asymmetry of the electronic dis- 
persion manifest itself in the WL behavior. 



2.2 Interference of electronic waves in monolayer graphene 



The WL correction to conductivity in disordered conductors is a result of the constructive 
interference of electrons propagating around closed loops in opposite directions [11] as sketched 
in figure [3]^b) . Such interference is constructive in metals and semiconductors with negligibly 
weak spin-orbit coupling, since electrons acquire exactly the same phase when travelling along 
two time-reversed trajectories. 

WL is usually described [11] in terms of the particle-particle correlation function, Cooperon. 
Following the example of Cooperons for a spin |, we classify Cooperons as singlets and triplets 
in terms of 'isospin' (AB lattice space) and 'pseudospin' (inter/intra- valley) indices (see Sec- 
tion [2TTJ) . In fact, with regards to the isospin (sublattice) composition of Cooperons in a disor- 
dered monolayer, only singlet modes are relevant. This is because a correlator describing two 
plane waves, ^k + , p and $k_,- p Eqs. (|3I4I) . propagating in opposite directions along a ballistic 
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segment of a closed trajectory as in figure E^b) has the following form: 

%,p%',-p ~ I T)k, p | l)K',-p - | 4>K, P | T)k',- p - T)k )P | T)k',- p + e lv \ I>k, p | I)k',- p . 

It contains only sublattice-singlet terms (the first two terms) because triplet terms (the last 
two terms) disappear after averaging over the direction of momentum, p = (p cos (p,p sinyj), so 
that (e ±l¥ % = 0. In fact, our diagrammatic calculation described in Section 12.41 shows that 
the interference correction to the conductivity of graphene is determined by the interplay of 
four isospin singlet modes: one pseudospin singlet and three pseudospin triplets. Of these, two 
of the pseudospin triplet modes are intravalley Cooperons while the remaining triplet and the 
singlet are intervalley Cooperons. 

In the WL picture for a diffusive electron in a metal, two phases $i and $2 acquired while 
propagating along paths "1" and "2" [see figure[3l(b)] are exactly equal, so that the interference 
of such paths is constructive and, as a result, enhances backscattering leading to WL [11]. In 
monolayer graphene the Berry phase 7r characteristic for quasi-particles described by the first 
term of Hi, determines the phase difference S = fii — 1?2 = kN (where N is the winding number 
of a trajectory) [10115] , and one would expect weak anti-localization behavior. However, the 
asymmetry of the electron dispersion due to hi w , leading to warping of the Fermi line around 
each valley as in figured]^), deviates S from ttN. Indeed, any closed trajectory is a combination 
of ballistic intervals, figure [3jb). Each interval, characterized by the momenta ±pj (for the two 
directions) and by its duration tj, contributes to the phase difference Sj = [e(pj) — e(—pj)]tj = 
hi w (jpj)tj. Since Sj are random uncorrelated, the mean square of S = ^ Sj can be estimated as 
(S 2 ) ~ ((tjhi w (pj)) 2 )t/T tr , where t is the duration of the path and Tt r is the transport mean 
free time. Warping thus determines the relaxation rate, 

r" 1 ~ (Trh 2 lw (p)) v , (5) 

which suppresses the two intravalley Cooperons, and, thus, weak anti-localization in the case 
when electrons seldom change their valley state. The two intervalley Cooperons are not affected 
by trigonal warping due to time-reversal symmetry of the system which requires e(K±,p) = 
e(K T , — p), figure 12b). These two Cooperons cancel each other in the case of weak intervalley 
scattering, thus giving Sg ~ 0. However, intervalley scattering, with a rate larger than 
the decoherence rate r^ 1 , breaks the exact cancellation of the two intervalley Cooperons and 
partially restores weak localization. 

2.3 Matrix parameterization, valley symmetry and the model of disorder 

To describe the valley symmetry of monolayer graphene and parameterize all possible types of 
disorder, we introduce two sets of 4x4 Hermitian matrices £ = (£ x , S y , S z ) with [S Sl , 2J S2 ] — 
2ie SlS2S3 S S3 , and 'pseudospin' A = {A x ,A y ,A z ) with [A h ,A h } = 2ie hhh defined as 

s x = n z ® a x , s y = ii z ® o y , E z = n Q ® a z , (6) 

A x = n x <g> cr 2 , A y =n v ® a z , A z = n z ® a . (7) 

The operators S and A form two mutually independent algebras equivalent to the algebra of 
Pauli matrices (in Eqs. (|6IT|) e SlS2S3 is the antisymmetric tensor and [i? S) yl;] = 0) thus they 
determine two commuting subgroups of the group U4 of unitary transformations |31j of a 4- 
component <P: an 'isospin' sublattice group SV2 = {e wnS } and a 'pseudospin' valley group 
SXJ2 = {e zbnA }. Also, S and A change sign under the inversion of time, whereas products 
S s Ai are invariant with respect to the t — ► — t transformation and can be used as a basis for 
a quantitative phenomenological description of non-magnetic static disorder [29132] . Table 1 
is a summary of the discrete symmetries of the operators S and A and their products S s Ai. 
Time reversal T of an operator W is described by (II X ® a x )W*{II x ® a x ). The rotation by 7r/3 
about the perpendicular z axis is described by Cq = TI X ® exp[(— 27ri/3)a z ]. Reflection in the 
x-z plane is R x = TTq <J X . 
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Fig. 3. (a) angular dependence w(ip) ~ cos 2 (y/2) of the scattering probability off a short range 
potential in monolayer graphene, (b) a pair of closed paths which contribute to weak localization, (c) 
angular dependence w(tp) ~ cos 2 (</?) of the scattering probability off a short range potential in bilayer 
graphene. 



The operators X and A help us to represent the electron Hamiltonian in weakly disordered 
graphene as 

Hi = v Sp + hi w + iu(r) + ^ £ s A/u S)i (r), (8) 

s,Z— x,y,x 

where /ii w = -fiS x ( J}p)A z E x ( Sp)^. 

The Dirac-type part v Sp of Hi in equation (jSj) and potential disorder Iu(r) (where I is a 
4x4 unit matrix and (u(r)u(r')) = u 2 <5 (r — r')) do not contain valley operators A;, thus, 
they remain invariant with respect to the pseudospin transformations from valley group SU2 • 
Below, we assume that the isospin/pseudospin-conserving disorder due to charges lying in a 
substrate at distances from the graphene sheet shorter or comparable to the electron wavelength 
H/pf dominates the elastic scattering rate, r _1 w Tq 1 — nyip/h, where 7 = pf/C^itTi v) 
is the density of states of quasiparticles per spin in one valley. All other types of disorder 
which originate from atomically sharp defects |29|32j and break the SU2 pseudospin symmetry 
of the system are included in a random matrix £ s Aiu s j(r). In particular, u z . z (r) describes 
disorder due to different on-site energies on the A and B sublattices, u x ( y y z (r) plays the role 
of a valley-antisymmetric vector potential of a geometrical nature, and u SiX m(t) take into 
account inter-valley scattering. For simplicity, we assume that different types of disorder are 
uncorrelated, (u s ,i(t)u s >,i' (r')) = u^iS S s'Sii'S(r — r') and, on average, isotropic in the x — y plane, 
u xi — u yi = u \n U L = u % = u s_l- We parametrize them by scattering rates T~ t — nju^/h. 
Also, the warping term, hiv, lifts the pseudospin symmetry SU2 , though it remains invariant 
under pseudospin rotations around the z-axis. 

To characterize Cooperons in monolayer graphene, we use a Cooperon matrix C^j^J^, where 
subscripts describe the isospin state of incoming a(3 and outgoing a' [3' pairs of electrons and 
superscripts describe the pseudospin state of incoming £/i and outgoing £'// pairs. Following 
the example of Cooperons for a spin i we classify Cooperons as singlets and triplets in terms 
of isospin and pseudospin indices C^? 1 ^ 2 . For example, M — is a 'pseudospin-singlet', M = 
x, y, z are three 'pseudospin-triplet' components; S = is a 'isospin-singlet' and S = x,y, z are 
'pseudospin-triplet' components. It is convenient to use pseudospin as a quantum number to 
classify the Cooperons in graphene because of the hidden SU^ symmetry of the dominant part 
of the free-electron and disorder Hamiltonian. 

2.4 Diagrammatic calculation of the weak localization correction in monolayer graphene 

To describe the quantum transport of 2D electrons in graphene we evaluate the disorder- 
averaged one-particle Green's functions, vertex corrections, calculate the Drude conductivity 
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Table 2. Matrices E a and Ai provide us with representations of the crystalline symmetry group, 
which is constructed of 3 generators, -—rotation, Ce, mirror reflection with respect to Ox axis, R x , and 
translation along Ox by lattice constant, a. Operation T stands for the time reversal, t — > —t. Here, E„ 
and Ai are grouped into bases forming irreducible representations which can be 1, 2 and 4 dimensional. 
The transformation matrixes Uji, gfa = Uji(g)fa< (g stands for a symmetry operation), are given 
explicitly for each of such bases fa. 
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and transport time, classify Cooperon modes and derive equations for those which are gapless 
in the limit of purely potential disorder. In Section 12.51 we analyse 'Hikami boxes' |ll|12j for 
the weak localization diagrams paying attention to a peculiar form of the current operator 
for Dirac electrons and evalute the interference correction to conductivity leading to the WL 
magnetoresistance. In these calculations, we treat trigonal warping /i lw in the free-electron 
Hamiltonian Eqs. Q2I8|) perturbatively, assume that potential disorder lu{r) dominates in the 
clastic scattering rate, t -1 sw t^ 1 = tt^u 2 /ft, and take into account all other types of disorder 
when we determine the relaxation spectra of low-gap Cooperons. 

Using the standard methods of the diagrammatic technique for disordered systems |11|12] 
and assuming that p-pvr ^> ft, we obtain the disorder averaged single particle Green's function, 

G R/A (P, e) = y A+ lfj , e R/A = e ± ^Hr \ (9) 

e R/A U P 

Note that, for the Dirac- type particles described in equation ((5J), the current operator is a 
momentum-independent matrix vector, v = t>S. As a result, the current vertex Vj ( j — x, y), 
which appears as a block in figure [DJa) describing the Drude conductivity, 

9ji = ^ / -|^Tr {vjG R (p, e) VjG A (p, e)} = ie^D; D = v 2 r = \v 2 r tll (10) 

is renormalised by vertex corrections [55] m figure [31(b) : v = 2v = 2vT,. Here 'Tr' stands for 
the trace over the AB and valley indices. Using the Einstein relation in equation (fTT)]) , we see 
that due to the anisotropy of scattering [i.e., lack of backscattering from an individual Coulomb 
centre as in figure [3^ a)] the transport time in graphene is twice the scattering time, r tr = 2tq. 
Note that in equation (|10|) spin degeneracy has been taken into account. 

The Cooperon C^'^J 1 ^, obeys the Bethe-Salpeter equation represented diagrammatically in 
figure fflje). The shaded blocks in figure 0Jc) are infinite series of ladder diagrams, while the 
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Fig. 4. (a) Diagram for the Drude conductivity with (b) the vertex correction, (c) Bethe-Salpeter 
equation for the Cooperon propagator with valley indices // and AB lattice indices af3a'/3' . (d) 
Bare 'Hikami box' relating the conductivity correction to the Cooperon propagator with (e) and (f) 
dressed 'Hikami boxes'. Solid lines represent disorder averaged G fl//A , dashed lines represent disorder. 



dashed lines represent the correlator of the disorder in equation ijHJ). Wc classify Cooperons in 
graphene as iso- and pseudospin singlets and triplets, as was mentioned above, with the help 
of the following relation, 

c i\% = \ E E (^.xA/^J^cl^C^^A)^- (n) 

a,/3,a',0', /u', 

Such a classification of modes is permitted by the commutation of the iso- and pseudospin 
operators X and A in Eqs. (|6I7I11[) . [27 s ,.Aj] = 0. To select the isospin singlet (s = 0) and 
triplet (s = x, y, z) Cooperon components (scalar and vector representation of the sublattice 
group SUf = {e ianS }), we project the incoming and outgoing Cooperon indices onto matrices 
E y £ Sl a,Tid E S2 S y , respectively. The pseudospin singlet (I = 0) and triplet (I = x,y,z) Cooper- 
ons (scalar and vector representation of the valley group SU^ 1 = {e lbn A }) are determined by the 
projection of C^'^/^, onto matrices A y Ai x (Ai 2 A y ) and are accounted for by superscript indices 

For 'diagonal' disorder Iu(r), the Bethe-Salpeter equation, figure [4|c) takes the form 
C L 6 \%(ci) =t S 1 ^5 s1S2 (12) 
+ I^E C ^( f l) / -^^{ZsZyAiAy [G^ Ru)+e f A y A k S y S Sl G^_ p>e }. 

s,l V ' 

It leads to a series of coupled equations for the Cooperon modes C l s l s = C l s . It turns out that 
for potential disorder Iu(r) isospin-singlet modes C l are gapless in all (singlet and triplet) 
pseudospin channels, whereas triplet modes C l x and C l y have relaxation gaps 7^ = F y = ^Tq 1 

and C l z have gaps F l z = Tq 1 . When obtaining the diffusion equations for the Cooperons using 
the gradient expansion of the Bethe-Salpeter equation we take into account its matrix structure. 
We find that isospin-singlets C l are coupled to the triplets C l x and C l y in linear order in the small 
wavevector q, so that the derivation of the diffusion operator for the isospin-singlet components 
would be incorrect if coupling to the gapful modes were neglected. The matrix equation for each 
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After the isospin-triplet modes are eliminated, the diffusion operator for each of the four 
gapless/low-gap modes C l becomes Dq 2 — iu; + Jq, where D = \v 2 T tT — v 2 t . 

Symmetry-breaking perturbations lead to relaxation gaps _Tq in the otherwise gapless pseudo- 
spin-triplet components of the isospin-singlet Cooperon Cq . All scattering mechanisms described 
in equation ([5]) should be included in the corresponding disorder correlator (dashed line) on 
the r.h.s. of the Bethe-Salpeter equation and in the scattering rate in the disorder-averaged 

QR/A 

as t 1 — * r 1 = t 1 + J2si T si 1 - This opens relaxation gaps in all pseudospin-triplet 
modes, Cq,Cq,C§, though does not generate a relaxation of the pscudospin-singlet Cq which 
is protected by particle conservation. 

The trigonal warping term hi w in the free electron Hamiltonian equation ((21) breaks the 
p — > p symmetry of the Fermi lines within each valley [33]. It has been noticed [34] that 
the deformation of a Fermi line of 2D electrons in GaAs/AlGaAs heterostructures in a strong 
in-plane magnetic field suppresses Cooperons as soon as the deformation violates p — > p 
symmetry. As h\ w has a similar effect, it enhances the relaxation rate of the pseudospin-triplet 
intravalley components Cq and Cq by 



2r (eV^ 2 )' 



(14) 



The estimated warping-induced relaxation time is rather short for all electron densities in the 
samples studied in [17], r w /r tr ~ 5 — 30, t w < t v , which excludes any WAL determined by 
intravalley Cooperon components. However, since warping has an opposite effect on different 
valleys, it does not lead to relaxation of the pseudospin-singlet Cq or the intervalley component 
of the pseudospin triplet, Cq. 

Altogether, the relaxation of modes Cq can be described by the following combinations of 
rates: 

r ° = 0, r* = 2rr\ rg = r% = r- 1 ee r w 1 + r" 1 + rr\ (15) 
where is the intervalley scattering rate (here we use the x — y plane isotropy of disorder, 



1 ~ t. i 1 and r 



I ~ 'yl — 'JJ )' 

rr 1 = 4rji + 2r~l, and r" 1 = 4rjJ + 2t zz \ 



(16) 



After we include dephasing due to an external magnetic field, B = rotA and inelastic decoher- 
ence, r" 1 , the equations for C read 



D(zV + HA)' +n + t- 1 ~ z W ]C (r, r') = S (r - r') 



(17) 



2.5 Weak localization magnetoresistance in monolayer graphene 

Due to the momentum- independent form of the current operator v = 2wS, the WL correction 
to conductivity Sg includes two additional diagrams, figure |4je) and (f ) besides the standard 
diagram shown in figure|3Jd). Each of the diagrams in figure^e) and (f) produces a contribution 
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equal to (—j) of that in figure B|d). This partial cancellation, together with a factor of four 
from the vertex corrections and a factor of two from spin degeneracy leads to 



5g 



2e 2 D 



~<1 



(CS + Ci> + C^-C° ) 



(18) 



Using equation (|18p . we find the B = temperature dependent correction, Sp/p = —5g/g, 
to the graphene sheet resistance. Taking into account the double spin degeneracy of carriers we 
present 



S P (0) 
P 2 



e 

nh 



ln(l + 2^) 



2 In 



T y /r tr 

1 + ^ 



and evaluate magnetoresistance, p{B) — p(0) = Ap(B), 
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e 2 p 2 
irh 
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B 



B v + £>* 
he 



(19) 



(20) 



Here, ip is the digamma function and the decoherence (taken into account by the rate r„ ) 
determines the curvature of the magnetoresistance at B < B v = Tic/ADer^. 

The last term in equation (fT5)) . Cq is the only true gapless Cooperon mode which determines 
the dominance of the WL sign in the quantum correction to the conductivity in graphene 
with a long phase coherence time, t v > Ti. The two curves sketched in figure [T] illustrate the 
corresponding MR in two limits: B* <C Bi (r* » Ti) and i?* ^> Bi (t* <C t ; ). In both cases, 
the low-field MR {B < B{) is negative. If B» < the MR changes sign: Ap(B) < at 
B < Bi = hc/ADeTi and Ap(B) > at higher fields. For B» > B h the MR is distinctly 
of a WL type, with almost no sign of WAL. Such behavior is expected in graphene tightly 
coupled to the insulating substrate (which generates atomically sharp scatterers). In a sheet 
loosely attached to a substrate (or suspended), the intervalley scattering time may be longer 
than the decoherence time, t\ > t v > t w (Bi < B v < BA. Hence Cq is effectively gapless, 
whereas trigonal warping suppresses the modes Cq and Cq. In this case the contribution from 
Cq cancels Cq, and the MR would display neither WL nor WAL behavior: Ap(B) = 0. 



3 Weak localisation magnetoresistance in disordered bilayer graphene 

3.1 Low energy Hamiltonian of bilayer graphene 

Bilayer graphene consists of two coupled monolayers. Its unit cell contains four inequivalent 
sites, A,B,A and B (A, B and A, B lie in the bottom and top layer, respectively) arranged 
according to Bcrnal stacking |35|5j : sites B of the honeycomb lattice in the bottom layer lie 
exactly below A of the top layer, figure \5\ The Brillouin zone of the bilayer, similarly to the 
one in monolayer, has two inequivalent degeneracy points K + and K_ which determine two 
valleys centered around e = in the electron spectrum [33] . Near the center of each valley 
the electron spectrum consists of four branches. Two branches describing states on sublattices 
A and B are split from energy e = by about ±71, the interlayer coupling, whereas two 
low-energy branches are formed by states based upon sublattices A and B. The latter can be 
described [5] using the Hamiltonian, which acts in the space of four-component wave functions 
<P = [</>k + ,a,0k + sj^k si < / ) k_,.a], where <j)^ t u is an electron amplitude on the sublattice a = 
A, B and in the valley £ = K + , K_. 

H 2 L = -g— [{Pi - Py) O x + 2p x p y a y ] + h 2w + V disorder , (21) 

h 2 w = v^n z (p x a x - Py(T y ) . 
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Fig. 5. Schematic of the bilayer lattice (bonds in the bottom layer A, B are indicated by solid lines 
and in the top layer A, B by dashed lines) containing four sites in the unit cell: A (white circles), B 
(grey), AB dimer (black). 



Here, cr x<ytZ and H Xl y.z are Pauli matrices acting in sublattice and valley space, respectively. 

The first term in equation (|21[) is the leading contribution in the nearest neighbors approx- 
imation of the tight binding model [5] . This approximation takes into account both intralayer 
hopping A <r-¥ B and A <-> B (that leads to the Dirac-type dispersion e = ±pv near the Fermi 
point K± in a monolayer) and the interlayer A <-> B hopping. This term yields the parabolic 
spectrum e = ±p 2 /2m with m = r yi/2v 2 which dominates in the intermediate energy range 
jji( v 3/ v ) 2 < £f < j7i- In this regime we can truncate the expansion of H(p) in powers of the 
momentum p neglecting terms of the order higher than quadratic. Electron waves characteristic 
for the first, quadratic, term of H2L have the form 

p ±ipx/R 

<?K, ±P = ±-^5- T)k,±p - e*| |)k, ±p ) , (22) 

where | T)k + ,± p = [1,0,0,0], | |) K+ ,± P = [0,1,0,0] and | T)k_,± p = [0,0,1,0], | T)k_,± p = 
[0, 0, 0, 1]. These are eigenstates of an operator ar\2 with a\\2 = — 1 for electrons in the conduc- 
tion band and crri2 = 1 for electrons in the valence band, where n2(p) = (cos(2y>), sin(2y>)) for 
p = (pcost/Jjpsint/?), which means that they are chiral, but with the degree of chirality different 
from the one found in monolayer (see Sec. 12. ip . Such electron waves are characterized by the 
Berry phase 27r, and the dependence w(9) ~ cos 2 9 of the scattering probability off a short-range 
potential on the scattering angle 9 — pp' is such that transport and scattering times in the 
bilayer coincide, although w(9) is anisotropic [see figure OUc)], and the Drude conductivity of a 
bilayer is g = 4e 2 nro/m (in contrast to monolayer graphene, see Sec. 12. ip . 

The second term in equation (|21[) . /i2w> originates from a weak direct A <-» B interlayer 
coupling. It leads to a Lifshitz transition in the shape of the Fermi line of the 2D electron gas 
which takes place when cf ~ £l = / v) 2 . In a bilayer with ep < cl, the interplay between 

the two terms in H^l determines the Fermi line in the form of four pockets [S] in each valley. In 
a bilayer with £f > cl, ^2w can be treated as a perturbation leading to a trigonal deformation of 
a single-connected Fermi line, thus manifesting the asymmetry of the electron dispersion inside 
each valley: e(K±,p) 7^ e(K±,— p). This asymmetry leads to the dephasing effect of electron 
trajectories similar to the one discussed in the case of monolayer, and is characterized by the 
scattering rate t" 1 equation (JSJ). 

The term Vdisorder in the equation (|21|) describes time-reversal-symmetric disorder. It is 
parameterized using t — ► — t symmetric 4x4 matrices acting in the sublattice/ valley space, 
which are listed in Table 3. 



Vdisorder ^ 

IIia s u s i{v) + II z a z u zz (r). 

s,l—0,x,y 



(23) 



12 



Will be inserted by the editor 



Table 3. Transformations of matrices of the form IIi<j 3 ,s,1 = 0,x,y,z, under crystalline symmetry 
group generators and time reversal operation. In bilayer graphene rotations and reflection symme- 
try operators are multiplied by the operation of reflection with respect to the z — plane, which is 
equidistant with respect to two honeycomb lattice layers. Therefore symmetry group generators are 
^-rotation, CeR z , mirror reflection with respect to Ox axis, R X R Z , and translation along Ox by lat- 
tice constant, a. Operation T stands for the time reversal, t — > — t. Here matrices are grouped into 
bases forming irreducible representations of the symmetry group which can be 1, 2 and 4 dimensional. 
The transformation matrixes Uji, gfa = . Uji(g)<f>j (g stands for a symmetry operation), are given 
explicitly for each of such bases <f>i. 
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The sum in equation (|23[) contains valley and isospin conserving disorder potential Iu(r), with 
(w(r)u(r')) = u 2 <5(r — r') and Tq 1 — T7"/u 2 /h : j = which originates from charged impurities 
in the SiC>2 substrate and is assumed to be the dominant mechanism of scattering in the system. 
All other types of disorder which breaks valley and sublattice symmetries are assumed to be 
uncorrelated, (u 8 i(r)u s i r(r')) = u 2 s1 5 ss i5w8(y — r'). We characterize them using scattering rates 
T^ 1 = TT^/ulJh. Furthermore, the scattering is assumed to be isotropic in the x — y plane, so 
that u xl — u yl = u ± i, u sx — u sy = u s± . 

3.2 Interference of electronic waves in bilayer graphene 

To analyze the WL effect we introduce Cooperon matrix C^ a t^, where subscripts describe the 
sublattice state of incoming a(3 and outgoing a 1 (3 1 pairs of electrons and superscripts describe 
the valley state of incoming £/i and outgoing pairs. Note that in contrast to monolayer we do 
not rewrite the bilayer Hamiltonian in terms of E and A matrices. We parametrize Cooperons 
as Cg 1 ^ 2 by M±,M2 "valley" and S\,S2 "sublattice" singlet and triplet states in a similar 
way to monolayer isospin and pseudospin states. The sublattice composition of Cooperons is 
determined by the correlator of plane waves propagating ballistically in opposite directions, 

*k,p*k',-p~| T)k, p I i>K',- P + | !>K,p| T)K',- P -e 2 ^| T)k, p | T)K',-p-e- 2 ^| l)K,p| DK'.-p. 

It is seen from the above expression that after averaging over the momentum direction the terms 
corresponding to C** y cx (| T)k, p | T)k',- p ±| 4)k, p | 4)k',- p ) disappear, since p = (p cos if, p simp) 
so that (e ±2l(p ) v = 0, whereas terms correponding to the sublattice symmetric Cooperons, 
Cf cx (| T)k, p | I)k',-p + I I)k, p | T)k',- p ) remain non-zero. 
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The dephasing effect of trigonal warping in bilayer is similar to monolayer, although it is 
caused by a different mechanism, its magnitude is estimated by equation ([5]). Dephasing due to 
warping suppresses the intravalley Cooperons Cf ,y leading to the absence of WL magnetore- 
sistance in the case of weak intervalley scattering, t\ 3> t v . In the case of strong intervalley 
scattering, t\ <C t v , WL is partially restored, thus, we predict the WL behavior of bilayer 
graphene with strong trigonal warping of Fermi line in each valley to be described by equa- 
tion ©. 



3.3 Diagrammatic calculation of the weak localization correction in bilayer graphene 

We derive the disorder averaged Green function for the bilayer Hamiltonian equation (j2Tj) : 



G«' A (p,e)= eR/A - e > <m ; (p ) 



e 2 _ £ 2 ' ( 24 ) 

where €r/a = e± ^ihr" 1 and r _1 = Tq 1 + t^ 1 + r" 1 w Tq . Here we introduced the following 
notations for the scattering rates r^ 1 = nju^/h, where t~ x — T~ y x = t^ 1 and r^ 1 = r^ 1 = tJj 

can be combined into the intervalley scattering rate Tj -1 = 4t7^ + 2t z ^ l 1 and the intravalley rate 
r^ 1 = r^ 1 both of which lead to an additional suppression of intravalley modes. Intervalley 
scattering also leads to the relaxation of C° although it does not affect the valley-symmetric 
mode Cf . Together, all the scattering mechanisms limit the transport time r" 1 = Tq +Yl s l t~, . 

Due to quadratic spectrum of quasiparticles in bilayer graphene the velocity operator, v x = 
— (jPxQx + Py& y )/m, v y — (p y a x — p x <jy)/m, is momentum dependent, and thus the current 
vertices in the conductivity diagram figure SJc) are not renormalized by impurity scattering 
accounted for by the diagram series figure 2Kb). As a result, the Drude conductivity is described 
by g = 4e 2 z/D, where D — \v 2 f tq and r tr = To. 

We parametrize the Cooperons utilizing the expression, 

c s&* = \ E E (wsXp in y n Ki r c^ : % (a S2 a y ) , a , (n K2 n v r'z' . 

The Bethe-Salpeter equation for Cooperons in bilayer reads, 

^SiS 3 " fa) = T ^Mi,M2^Si,S 2 

+ i^E C 5Sf 2 fa) J j^Tr{n M a s n y a y [G* hu)+£ Y n y a y n Ml a Sl Gh_ Vie }.(25) 

We find that Cgg}' 1 ' = 5 MM 5ss'C<f and that sublattice-singlet has a relaxation gap 
-To 1 = Tq , sublattice-triplets C£f , have gaps r^ 1 = r* f = \t^ x , whereas symmetric 

sublattice-triplet Cooperon C^ 1 is gapless. Due to warping of the Fermi line induced by hi w in 
the free-electron Hamiltonian ([2"Tjl . the intravalley Cooperons , C| are suppressed, even in 
a bilayer with purely potential disorder. Warping opens a gap, r^ 1 in the relaxation spectrum 
of these 'valley-triplet' Cooperon components: 

t-1 - J 2^ T ( Tr ^w(p))v = ttolIV- 1 , nn L l 2 < 1 . . 

w ~ 1 t-\ nn L l 2 > 1 ' 1 ° J 

where til is the density of electrons at which Lifshitz transition occurs, I and r are mean free 
path and transport time in the system respectively. We estimate that for the recently studied 
bilayers [9] with n e = 2.5 x 10 12 cm~ 2 , r w ~ r and the mean free path I ~ 0.1/xm. A similar 
situation occurs in bilayer structures studied by R. Gorbachev et al. [18]. Also, a short-range 
symmetry-breaking disorder w.y generating intervalley scattering leads to the relaxation of C_j, 
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although it does not affect the valley-symmetric mode Cf. Thus we find that the low-gap modes 
C^ 1 obey the diffusion equation, 



r + r-i + z^v + iA) 2 



C(r,r') = 5(r-r'), 



r z 



o, rl = 2rr\ 



(27) 



px{y) -1 _ 

z ' * — 



2r; 



where we included dephasing due to an external magnetic field, B = rotA, temperature- 
dependent inelastic decoherence, r~ (T), and all of the above mentioned relaxation mecha- 
nisms. 



3.4 Weak localization magnetoresistance in bilayer graphene 

The interference correction to the conductivity in a bilayer can be expressed in terms of C (r, r), 
the solutions of the above Cooperon equations taken at coinciding coordinates: 

5 9 = ^ + C° z - CI - C*] . (28) 

For completeness, in equation (f2"8"|) we have retained the intravalley Cooperons ' v , though they 
are strongly suppressed by trigonal warping. Following their suppression, the WL correction is 
determined by the intervalley modes C® and Cf but, in the absence of intervalley scattering, the 
contributions of C° and Cf are equal in magnitude, so that they cancel. Intervalley scattering 
due to atomically sharp scatterers breaks this exact cancellation and partially restores the WL 
effect. Equations (|28|27p yield the zero field WL correction to the resistivity and the WL MR, 



, In 1 + 2-* +<*,,, (29) 
p nil \ Tj 



Ap(B) _ e 2 p 
p irh 



' , B s , B 
F(—)-F(- 



B v B v + 2Bi ' 



6(B), 



where B Vtl = Uc/(4DeT Vi i). Equation (f2"9"| gives a complete description of the crossover be- 
tween two extreme regimes mentioned at the beginning [25] . It also includes small contri- 
butions of the suppressed intravalley Cooperons, <5 = [2e 2 p/(wh)] ln(r v T*/[r(r* + r v )]) and 
5(B) = -[2e 2 p/(TTh)]F[B/(B ip + where r" 1 = r" 1 + 2T" 1 + rf 1 an d B, = fic/(4DeT„). 
This permits us to account for a possible difference between the warping time t w and the trans- 
port time r. According to equation (f2"9"]l WL MR in bilayer graphene sheet disappears as soon 
as T; exceeds t v , whereas in structures with t v > t\, the result equation (|29p predicts the WL 
behaviour, as observed in [18]. Such WL MR is saturated at a magnetic field determined by the 
intervalley scattering time, instead of the transport time as in usual conductors, which provides 
the possibility to measure r; directly. 



4 Conclusions and the effect of edges in a disordered nanoribbon 

We have shown that p — > — p asymmetry of the electron dispersion in each valley of graphene 
leads to unusual (for conventional disordered conductors) behavior of interference effects in 
electronic transport. Without intervalley scattering, trigonal warping of the electron disper- 
sion near the center of each valley destroys the manifestation of chirality in the localization 
properties, resulting in a suppression of weak anti- localization in monolayer graphene and of 
weak localization in a bilayer. Intervalley scattering tends to restore weak localization, and this 
behavior is universal for monolayer and bilayer graphene, despite the fact that electrons in 
these two materials have different chiralities and can be attributed different Berry phases: it in 
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monolayers, 2n in bilayers |3|5j . This suggests that a suppressed weak localization magnetore- 
sistance and its sensitivity to intervalley scattering are specific to all ultrathin graphitic films 
independently of their morphology [19] and are determined by the lower (trigonal) symmetry 
group of the wavevector K in the corner of the hexagonal Brillouin zone of a honeycomb lattice 
crystal. 

The influence of intervalley scattering on the WL behavior determines a typically negative 
(WL) MR in graphene nanoribbons. Indeed, in a narrow ribbon of graphene, monolayer or 
bilayer, with the transverse diffusion time L^/D <C t\,t*,t v , the sample edges determine 
strong intervalley scattering rate [36] . Thus, when solving Cooperon equations in a wire, we 
estimate ig ~ tt 2 D/L\ for the pseudospin triplet, whereas the singlet Cq remains gapless. This 
yields negative MR persistent over the field range B < 2ttB±, where B± = hc/eL 2 ^: 



Ap wirc (B) 2e% 



1 



1 + \B 2 /B ip B A 



(30) 



The results of Eqs. (|19I20I28I29|) . and (J3D1) give a complete description of the WL effect in 
graphene and describe how the WL magnctoresistance reflects the degree of valley symmetry 
breaking in it. 

This project has been funded by Lancaster-EPSRC Portfolio Partnership grant EP/C511743 
and was completed during the MPI PKS Seminar "Dynamics and Relaxation in Complex 
Quantum and Classical Systems and Nanostructures." 
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